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ABSTRACT 

Existing  rotadonally  invariant  electroelastic  equations  are  reduced  to  the  case  of  large 
electric  fields  and  small  strain.  These  latter  equations  are  specialized  to  the  case  of  thin  plates 
with  completely  electroded  major  surfaces,  and  it  is  shown  that  in  this  case  the  charge  equation 
of  electrostatics  is  satisfied  trivially  to  lowest  order.  It  is  further  shown  that  for  the  thin 
stress-free  polarized  ferroelectric  ceramic  plate  subject  to  large  electric  fields,  the  resulting 
equations  readily  account  for  experimental  data  in  existence  in  the  literature 
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1.  Introduction 

In  the  Held  of  active  vibration  control  thin  fully  electroded  polarized  ceramic  plates  are 

used  as  actuators  and  sensors.  The  equations  of  linear  piezoelectricity  are  employed  in  the 

analytical  description  of  the  polarized  ceramic  even  though  large  electric  driving  fields  are 

employed  during  actuation,  for  which  the  material  behavior  is  strongly  nonlinear.  To  account 

for  this  inherent  nonlinearity  it  has  been  suggested^  that  the  coefficients  of  the  linear  theory  be 

made  field  dependent,  which  is  inconsistent  from  a  theoretical  standpoint.  This  is  the  current 

state  of  affairs  in  spite  of  the  fact  that  a  rotationally  invariant  general  nonlinear  system  of 

electroelastic  equations,  which  can  readily  account  for  the  observed  behavior,  has  been  in 

2  3 

existence  in  the  literature  ’  for  many  years.  In  addition,  all  the  work  on  piezoelectric 
actuators  uses  only  the  linear  piezoelectric  constitutive  equations  and  does  not  discuss  the 
extent  to  which  the  stress  equations  of  motion  or  charge  equation  of  electrostatics  are  or  are 
not  satisfied. 

In  this  paper  the  general  rotationally  invariant  nonlinear  electroelastic  equations  are 
presented  and  specialized  to  the  case  of  large  fields  and  small  strain.  These  latter 
three-dimensional  equations  are  then  reduced  to  the  case  of  anisotropic  plane  stress  for  the  thin 
electroelastic  plate.  The  equations  for  the  thin  plate  are  further  specialized  to  the  case  when 
the  major  surfaces  are  completely  covered  with  electrodes,  and  it  is  shown  that  the  charge 
equation  of  electrostatics  is  satisfied  trivially  to  lowest  order  for  this  configuration.  The 
reduced  form  of  the  constitutive  equations  for  the  case  of  anisotropic  plane  stress  when  the 
large  surfaces  of  the  thin  plate  are  completely  covered  with  electrodes  is  obtained  from  the 
more  general  equations.  These  latter  equations  may  be  used  in  the  extensional  equations  of 
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motion  to  obtain  a  determinate  system  of  equations.  The  reduced  constitutive  equations  are 
specialized  to  the  case  of  the  polarized  ceramic  poled  normal  to  the  major  surfaces  and  it  is 
shown  that  the  resulting  equations  readily  account  for  measurements  by  Crawley  and 
Anderson^  of  the  stress- free  PZT  ceramic  plate  subject  to  large  electric  fields. 

2.  Electroelastic  Equations 

The  stress  equations  of  motion  and  charge  equation  of  electrostatics  for  material  points 
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of  an  electroelastic  solid  may  be  written  in  the  respective  forms 


^4. = »• 

where  we  have  employed  the  conventions  that  capital  and  lower  case  indices,  respectively, 
refer  to  the  reference  Cartesian  coordinates  and  present  Cartesian  coordinates  of  material 
points.  We  also  employ  the  conventions  that  a  comma  followed  by  a  capital  index  denotes 
partial  differentiation  with  respect  to  the  known  reference  coordinates,  i.e.,  the  independent 
variables  excluding  time,  a  dot  over  a  variable  denotes  partial  differentiation  with  respect  to 
time  and  repeated  tensor  indices  are  to  be  summed.  The  symbols  p°,  Uj^,  Ky,  and 
respectively,  denote  the  reference  mass  density,  the  mechanical  displacement,  the  total 
Piola-Kirchhoff  stress  tensor  (material  plus  Maxwell  electrostatic)  and  the  reference  electric 
displacement  vector,  which  is  related  to  the  ordinary  electric  displacement  vector  Dj,  by 

In  (2.3)  we  have  made  use  of  the  fact  that  the  motion  of  a  material  point  is  described  by  the 
relation 

y.  =  y.(XL.t),  (2.4) 

which  is  one-to-one  and  differentiable  as  often  as  required,  where  y^  and  Xj^,  respectively, 
denote  the  present  and  reference  position  of  material  points, 

J  »  d«  Xj 


(2.5) 


J 


and  a  comma  followed  by  a  lower  case  index  denotes  partial  differentiation  with  respect  to  the 
unknown  present  coordinates  of  material  points.  The  symbol  in  (2.1)  is  a  translation 
operator,  which  serves  to  translate  a  vector  from  the  present  to  the  reference  position  and 
vice-versa  and  is  required  for  notational  consistency  and  clarity  because  of  the  use  of  capital 
and  lower  case  indices,  respectively,  to  refer  to  the  reference  and  present  positions  of  material 
points.  Clearly,  the  mechanical  displacement  vector,  Uj^,  and  present  and  reference  positions 
of  material  points,  y^  and  Xj^,  are  related  by 


^i  ~  ^iM^^M 


(2.6) 


The  total  Piola-Kirchhoff  stress  tensor  and  reference  electric  displacement 
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vector  3i^  may  be  written  in  the  respective  forms 


•^Lj  =  «Lj  ^  My.  (2.7) 

^  i  * 

where  is  the  electric  permittivity  of  free  space,  MKS  units  are  being  employed  and 

^Lj  “  ~  ' 

The  symmetric  Maxwell  electrostatic  stress  tensor  material  strain  tensor  and 

rotationally  invariant  electric  variable  are  given  by 

"  2  '''l  ' 

The  electric  field  is  related  to  the  quasistatic  electric  scalar  potential  by 

E;  =  - 

and  ;i:  =  ;t(EKL’^L^  is  a  scalar  state  function,  which  for  purposes  of  tins  work  may  be  written 
in  the  form 


(2.11) 

(2.12) 

(2.13) 


~  I  ^ABCD^AB^CD  ‘  ®ABc''^A^BC  '  2  ^AB^A^B 
*  2  ^ABCD^A^B^CD  ’  5  ^ABC^A^B^C 

-  s^abcdeWaWb^cEde  ■  ^  Jabcd'^a^b^c'^D’  (2.14) 
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where  the  material  constants  c^gcD’  ®ABC’  |aB’  ^ABCD’  |aBC’  ^ABCDE’  JaBCD 

are  called  the  second  order  elastic,  piezoelectric,  second-order  electric  permeability, 
electrostrictive,  third-order  electric  permeability,  third  odd  electroelastic,  and  fourth-order 
electric  permeability,  respectively.  At  this  point  we  note,  from  (2.12)2  and  (2.13),  and  the 
chain  rule  of  differentiation,  that  the  rotationally  invariant  electric  vector  may  be  written 
in  the  simple  linear  form 

Wl  =  -4l-  »>5) 

The  system  of  equations  may  readily  be  reduced  to  four  equations  in  the  four  dependent 
variables  Uj^  and  (j)  by  straightforward  substitution.  The  four  equations  consist  of  (2.1)  and 
(2.2). 

The  boundary  conditions  across  material  surfaces  of  discontinuity  may  be  written  in  the 

form 

NL[KLjl=0,  luj^]=0,  (2.16) 

[(|»]=0,  (2.17) 

where  we  have  introduced  the  conventional  notation  [a  ]  for  at  -  aj,  N,  denotes  the  unit 
normal  directed  from  the  -  to  the  -h  side  of  the  reference  position  of  the  material  surface  of 
discontinuity,  and  3,  is  the  electric  surface  charge  density  per  unit  reference  area  that  can  exist 
at  a  surface  separating  an  insulator  from  a  conductor.  At  a  material  surface  separating  two 
insulators  3  vanishes.  The  boundary  conditions  may  readily  be  expressed  in  terms  of  the  four 
dependent  variables  and  ^  on  each  side  of  the  material  surface  of  discontinuity  by 

4 

straightforward  substitution  . 


systematic  treatment  of  this  and  the  next  section.  However,  we  note  that  in  the  last  section  of 
this  paper  we  present  the  extended  results  that  would  have  been  obtained  for  the  simple  case 
treated  there,  for  which  only  one  term  of  each  of  these  material  tensors  occurs. 

In  making  the  reduction  to  infinitesimal  strain  powers  of  Uj^  and  its  gradients  higher 
than  the  first  are  dropped  in  all  expressions,  and  even  the  linear  terms  themselves  are  dropped 
in  comparison  with  any  finite  quantity,  such  as  the  Kronecker  delta  or  1,  in  the  usual  manner. 
In  this  way  we  obtain 


(3.1) 

\m  ^LM  ^  2  ^“l,M  “M,L^- 

(3.2) 

Then,  introducing  (3.1)  and  (3.2)  in  (2.9)  and  (2.10)  and  employing  (2.11)  and  (2.15),  we 

obtain 

^Lj  ""  ^jMf^LMAB^AB  ®ALM*^,A  *  I  ^ABLM‘J*,A*^,B 

(3.3) 

*  ®LBC^BC  ■  ^AL^w\  2  ^ABL‘^‘,A^,B’ 

(3.4) 

^Lj  =  ^o^jM[‘^,L^,M  *  2  ’ 

(3.5) 

(3.6) 

The  substitution  of  (3.3)  and  (3.5)  into  (2.7)  enables  us  to  write 

^Lj  ^ 

(3.7) 

where 

\m  2LMAB^AB  ®ALM^,A  ’  2  ^ABLM'^’a'^.B’ 

(3.8) 

and 

^ABLM  ^  ^ABLM  ‘  ^^o[^^MB  ‘  2  ^AB^m]  ' 

(3.9) 
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The  further  substitution  of  (3.4)  and  (3.6)  into  (2.8)  yields 

""  ®LBC^BC  ■  ^LA***^  2  |ABL‘t’^‘t’,B’ 
where  the  tensor  of  dielectric  constants  is  given  by 

^A  ^  ^o^A  |lA- 

The  substitution  of  (3.7)  into  (2.1)  yields  the  stress  equations  of  motion  in  the  form 

and  we  still  have  the  charge  equation  of  electrostatics  (2.2),  which  we  rewrite  here  for 
completeness 


=  (3.13) 

At  this  point  we  have  the  stress  equations  of  motion  (3.12)  and  charge  equation  of 
electrostatics  (3.13),  both  referred  to  the  known  reference  coordinates.  These  are  the  four  field 
equations  for  the  electroelastic  continuum.  The  associated  constitutive  equations,  which  are 
linear  in  the  infinitesimal  strain  and  up  to  quadratic  in  the  electric  field  are  given  by  (3.8)  and 
(3.10).  The  linear  portion  of  these  constitutive  equations  is  identical  with  those  of  linear 
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piezoelectricity  ,  and  we  note  that  the  effective  electrostrictive  tensor  in  (3.8) 

contains  the  influence  of  the  Maxwell  electrostatic  stress  tensor  as  shown  by  (3.9)  with  (3.7) 
and  (3.5).  The  substitution  of  (3.8)  and  (3.10)  with  (3.2)  into  (3.12)  and  (3.13)  yields  four 
differential  equations  in  the  four  dependent  variables  Uj^  and  ({). 

4.  Extensional  Equations  for  the  Thin  Plate  with  Fullv  Electroded  M^or  Surfaces 

A  diagram  of  the  thin  plate  with  fully  electroded  major  surfaces  is  shown  in  Figure  1 
along  with  the  associated  coordinate  system.  Since  we  have  restricted  the  description  to 
infinitesimal  strain  in  Section  3,  we  dispense  with  the  no  longer  necessary  convention  that 
capital  indices  refer  to  the  reference  coordinates  and  lower  case  indices  refer  to  the  present 
coordinates  of  material  points,  and  use  lower  case  indices  exclusively  to  refer  to  the  reference 
coordinates  of  material  points,  to  which  the  entire  description  is  referred.  In  addition,  for 


/ 


convenience  and  because  it  is  conventional  we  introduce  the  usual  compressed  notation  for 
tensor  indices  for  stress  and  strain  in  the  constitutive  equations,  in  which  p,  q,  r,  s  take  the 
values  1,  2,  3,  4,  5,  6  as  ij  or  take  the  values  11,  22,  33,  23  or  32,  31  or  13,  12  or  12. 
Accordingly,  we  rewrite  the  constitutive  equations  (3.8)  and  (3.10)  in  the  compressed  notation 
thus 


where  we  have  taken  the  liberty  of  writing  for  and  for  since  the  distinction 
between  the  quantities  is  lost  when  the  assumption  of  infinitesimal  strain  is  employed. 

Since  the  tractions  on  the  electroded  major  surfaces  vanish  and  the  electroelastic 
plate  is  very  thin,  we  may  make  the  assumptions  of  plane-stress  and  take  =  0  throughout 
the  plate.  Similarly,  since  the  electric  fields  vanish  in  the  conducting  electrodes  and  the 
electroelastic  plate  is  thin,  we  may  take  E,,  =  0  throughout  the  plate,  where  we  have  introduced 

a 

the  convention  that  a,  b,  c,  d,  take  the  values  1  and  2,  but  not  3.  As  in  the  case  of  anisotropic 

plane  stress,  when  so  many  of  the  stress  components  vanish  it  is  advantageous  to  obtain  the 

strain-stress  relations.  To  this  end  we  operate  on  (4.1)  with  c^  =  s^  to  obtain 

S  =s  T  +d,  E.  +  i  E-E. , 
q  qp  p  kq  k  2'^kq  j  k’ 

where 


(4.3) 


^kp  ®kp^pq’  ^ktq  \^p^pq’ 

and  we  note  that  the  dj^^  are  a  well-known  set  of  previously  defined^  piezoelectric  constants, 
which  are  useful  for  plane-stress  tyjw  calculations,  whereas  the  plane-stress  type  nonlinear 


electroelastic  constants  are  new.  The  substitution  of  (4.3)  into  (4.2)  yields 

^Ik^k  Z  ^Tjk^j^k  ’ 


(4.5) 


where 
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(4.6) 


®  ^q^kq*  ^  ^ jk  "  ^/jk  ®  ^q^jkq’ 

T  6 

and  we  note  that  the  are  a  well-known  set  of  previously  defined”  dielectric  constants 

T 

while  the  nonlinear  electroelastic  susceptibility  constants  Xi^ 


In  the  present  notation  the  charge  equation  of  electrostatics  takes  the  form 


(4.7) 


Since  for  the  thin  plate  shown  in  Figure  1  the  D.  vary  slowly  with  respect  to  X  ,  we  first 

K  3 

rewrite  (4.7)  in  the  form 

°3,3*Da,a  =  »- 

and  then  note  that  D  is  of  smaller  order  than  D,  Hence,  to  lowest  order  we  have 

a, a 

D3  3  =  0,  (4.9) 

which  yields 

D3=D3(X^,t).  (4.10) 

Furthermore,  since 

E^  =  0,  (4.11) 

and  for  a  state  of  plane-stress  to  lowest  order  and  are  independent  of  X3,  we  have 


shown  that 


E3  =  -V/t,  (4.12) 

where  V  is  the  driving  voltage  across  the  electrodes  and  t  is  the  thickness  of  the  plate.  Thus 
we  have  shown  that  (4.11)  and  (4.12)  assure  that  (4.7)  is  satisfied  to  lowest  order  for  the 
nonlinear  constitutive  equations  (4.5)  or  (4.2).  Consequently,  with  (4.11)  and  (4.12),  Eq.(4.7) 
may  be  ignored  from  here  on  because  it  has  been  satisfied  to  lowest  order  for  the  nonlinear 
constitutive  equations  (4.2)  or  (4.5). 

Since  we  have  (4.11)  and  by  virtue  of  the  validity  of  plane  stress  for  the  thin  plate,  we 
also  have 


(4.13) 


the  constitutive  equations  (4.3)  and  (4.5)  may  be  written  in  the  reduced  form 

^qv'^v  ^3q^3  2  ^33q^3’ 

^  ^^3%  2  ^33 ^^3’ 


(4.14) 

(4.15) 


where  we  have  introduced  the  convention  that  v  takes  the  numbers  1,2,6  but  not  3,4,5.  For  the 
symmetry  of  the  polarized  ceramic  poled  in  the  X^-direction,  the  constitutive  equations  (4.14) 
and  (4.15)  when  written  out  in  full  take  the  reduced  form 

^1  ^  ^11^1  ^12^2  ^31^3  ^31^3’ 

^2  ^  ^12’^1  ^1 1’^2  ‘^31^3  2  hl^3' 

S3  =  Sj3Tj  +  SJ3T2  +  d33E3  +  j  ^33E^, 


S4  =  0,  -  0,  Sg  -  Dj  -  0,  D2  -  0, 

D3  =  d3j(Tj  +  T2)  +  e33E3  +  j 


(4.16) 


At  this  point  we  note  that  since  the  73^^^  have  been  taken  to  vanish  for  the  thin  plate, 

U3  must  be  negligible  on  account  of  (3.12),  which  is  always  the  case  for  the  low  frequency 
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extension  of  thin  plates  .  As  a  consequence,  from  (3.12)  the  stress  equations  of  motion  for 


extension  of  the  thin  plate  take  the  form 


Tab,a  = 

in  our  notation.  Since  for  the  low  frequency  extension  of  thin  plates  Eqs.(4.17)  remain  to  be 
satisfied,  we  must  solve  the  constitutive  equations  (4.16)  for  in  terms  of  S^j^  and  E3.  The 


result  is 


T  =cPs  +cPs  -e^E  -i^E^ 

-CiiOi +012^2  631^3  2*31^3’ 

^2  ^  ^12^1  ^11^2  '  ®31^3  ■  2  ^1^3’ 

^6  ^  ‘^66^6-  ^3  "  ®31^1  ®31^2  ^3^3  2  ^33^3’ 


(4.18) 


where  the  plate  constants  with  the  superscript  p  are  given  by 
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cP  =!y. 


P  _  12 


eP  =- 
®31  s 


*31 


11 


*12 


.  ^1  = 


*31 


+  s,^’  31  s 


11  ^2 


A  =  s2  -s^ 
A-Sii  5^2- 

£?  -fT  .  “31 

^33  "^33  s 


11  ^12 


p  ^  T  ;  '*31<^31 
^33  ^33  ^  s, ,  +  s 


‘11 


12 


(4.19; 


and  %vf  note  that  although  the  constants  e^j  and  l^^owr  the  expressions  for 

the  nonlinear  constants  and  substitution  of  (4.18)  with  (3.2)  into  (4.17) 

yields  tw'o  differentia]  equations  in  the  two  dependent  variables  u^.  When  the  major  surfaces 
of  the  polarized  ceramic  plate  are  completely  covered  with  conducting  electrodes,  the 
differential  equations  are  independent  of  the  voltage  V,  which  appears  only  in  certain  of  the 
boundary  or  continuity  conditions  at  the  edges. 


to  Large  Electric  Fields 
When  the  thin  ferroelectric  ceramic  plate  polarized  normal  to  the  major  surfaces,  which 
are  fully  electroded  and  subject  to  a  large  static  voltage,  is  stress-free,  the  stress  equations  of 
equilibrium  (4.17)  with  =  0  are  satisfied  trivially  and  the  constitutive  equations  for  the 
in-plane  strains,  (4.16)^  2>  relevant  form 


S]  -  S2  -  d3jE2  +  2 


(5.1) 


In  Section  3  we  ignored  the  material  tensors  ai'd  ^ABCD  systematic 


treatment  presented  in  Sections  3  and  4,  but  we  noted  that  we  would  present  the  extended 
result  for  the  relevant  constitutive  equations  that  would  have  been  obtained  for  the  simple  case 
treated  had  they  been  included.  For  the  stress-free  plate  the  constitutive  equations  that  would 
have  been  obtained  in  place  of  (5.1)  arc  given  by 

Sj  =  S2  =  d3iE3  +  2  ^31^3 [+  “6  >*331^3]  •  (5-2) 


1  1 


in  which  we  have  enclosed  the  term  that  was  not  obtained  in  the  systematic  treatment  in 
square  brackets. 

In  Figure  2  the  discrete  solid  boxes  denote  the  data  obtained  by  Crawley  and 

Anderson^  and  the  straight  dotted  line  tangent  to  the  data  near  the  origin  is  the  prediction  of 

linear  piezoelectricity  for  the  known  value  of  d^j  in  nnVV.  The  curved  solid  line  that  departs 

-4 

from  the  data  for  values  of  strain  larger  than  10  was  obtained  using  the  constitutive  equation 

given  in  (5.1)  with  the  known  value  of  d22  and  employing  a  least  squares  fit  to  the  measured 

data  for  values  of  microstrain  below  100  to  evaluate  .7949  pm  [V  .  The  curved  solid 

line  that  goes  through  all  the  data  points  was  obtained  using  the  constitutive  equation  given  in 

(5.2)  with  the  known  value  of  d^j  and  employing  a  least  squares  fit  to  all  the  measured  data  to 

evaluate  and  7332.  which  in  this  fit  have  the  values  .8054  pm  fW  and  7^32  =  ‘ 

3  3 

.7754  fin  /V  .  Thus  we  have  shown  that  the  general  nonlinear  electioelastic  description, 
which  existed  long  before  the  measurements  were  made,  can  readily  be  reduced  to  the  form 
that  describes  the  thin  electroelastic  polarized  ceramic  plate  subject  to  large  electric  fields  and 
accurately  accounts  for  the  measurements. 
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